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Abstract: We perform holographic renormalization for probe branes in AdS§ x S 5 . We 
show that for four known probe D-branes wrapping an AdS m x S n , the counterterms needed 
to render the action finite are identical to those for the free, massive scalar in AdS m plus 
counterterms for the renormalization of the volume of AdS m . The four cases we consider are 
the probe D7, two different probe D5’s and a probe D3. In the D7 case there are scheme- 
dependent finite counterterms that can be fixed by supersymmetry. 
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1. Introduction 

Holography equates a (d + l)-dimensional theory of gravity with a d-dimensional field the¬ 
ory. The AdS/CFT correspondence [1] is an explicit example of this, in which supergravity 
(SUGRA) formulated on the ten-dimensional bulk spacetime AdS$ x S 5 is dual to N = 4 
SU(N) super Yang-Mills (SYM) theory living on the four-dimensional boundary of the AdS§ 
factor. The precise statement of the duality equates the on-shell bulk action with the gener¬ 
ating functional of the CFT [2, 3]. In general, however, both of these quantities are infinite: 
the bulk action suffers IR divergences while UV divergences appear on the CFT side. The 
method of holographic renormalization [4, 5, 6, 7] was developed to remove these divergences 
in a consistent way so that the correspondence equates finite, physical quantities. 
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Holographic renormalization involves regulating the on-shell bulk action and then adding 
covariant counterterms to cancel the divergences that appear as the regulator is removed. 
This technique is very general, and should certainly apply to the subset of bulk solutions 
that describe supersymmetric (SUSY) probe brane solutions in AdS/CFT. The AdS/CFT 
duality is constructed from the near-horizon limit of a stack of N D3-branes with N —> oo. 
Introducing a finite number of additional branes orthogonal to the D3’s in some directions will 
produce probe branes embedded in the near-horizon geometry of the D3’s, which is AdS 5 x S 5 . 
The simplest known probe branes of this type preserve half the supersymmetry and introduce 
additional fields in the J\f = 4 SYM. Their action is the Dirac-Born-Infeld (DBI) action, whose 
contribution to the bulk action is simply the volume of the brane. Since the backreaction of 
the probe branes on the background geometry can be neglected this system is probably the 
simplest application of holographic renormalization in a string theory setting and serves as a 
nice illustration of the method. One aspect of holographic renormalization for submanifolds, 
the conformal anomaly, was previously studied in [ 8 ]. In this paper, we consider four known 
probes: a D7-brane [9], two D5-branes [10, 11], and a probe D3-brane [12]. One of the cases 
in non-supersymmteric [11] while the other three cases are half supersymmetric. The brane 
embedding in all of these cases is described by a single scalar field that, via the correspondence, 
will be dual to some operator in the field theory. 

Given a renormalized bulk action, renormalized correlators in the dual field theory can 
be computed straightforwardly. The goal of this paper is to compute renormalized one-point 
functions, or vacuum expectation values (vev’s), for the operators dual to the embedding 
scalars, using holographic renormalization. Especially for the D3/D7 system the situation is 
interesting since the vev is not just the subleading term in the asymptotic expansion of the 
scalar field. Our calculation of the vev clarifies the somewhat ad hoc procedure that has been 
used in the literature. 

This paper is organized as follows. In section 2, we review holographic renormalization, 
including the methods of adding covariant counterterms to render the bulk action finite and 
of computing one-point functions from this renormalized action. In section 3, we show that 
the counterterms for our DBI actions are identical to the counterterms for a free scalar in 
AdS. In section 4, for the D7, we use these counterterms to compute the one-point function 
for the embedding scalar’s dual operator. In this case, we find that a finite counterterm 
must be included if the renormalization scheme is to be supersymmetric. We do a similar 
computation of a one-point function in section 5 for the D5 probe. We also consider a second 
D5 embedding that is not considered in the general arguments of section 3 but which needs no 
new counterterms. In Section 6 we compute one- and two-point functions for the D3 probe. 


2. Review: Holographic Renormalization 

We review in this section the holographic renormalization of the volume of an AdS^+i space- 
time and of a free scalar in that spacetime. The former has been studied extensively in 
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[4, 5, 6 , 13], while the latter was explained in detail in [7], see [14] for a review of holographic 
renormalization. 

To begin, write the metric of asymptotically AdSd+i spacetime in the form of Fefferman 
and Graham [15] in units where the AdS radius is one (which we use throughout): 

dr 2 1 

ds 2 = G^dx^dx 1 ' = H— ^gij(x,r)dx l dx : ’ , ( 2 . 1 ) 

where r is the radial coordinate, fi, u = 0 ,..., d and i,j = 0,...,d — l and where 

g(x, r ) = s (0) + r 2 g [2) + ... + r d g {d) + h (d) log(r 2 ) + 0{r d+1 ). ( 2 . 2 ) 

The logarithmic term only appears for even d and only even powers of r appear up to order 
r d ~ l . Solving Einstein’s equation recursively for the g(,\ gives two useful identities 1 : Tr g^ 2 ) = 
with Rq the Ricci scalar built from g^, and Tr g^ = ^Tr (g 2 2 ^) [7]. 

For the renormalization of the volume of AdS d+ 1 , we restrict our discussion to d = 2,3,4, 
which will be relevant for the probe branes we consider. The volume is 2 

V AdS = J d d+1 xVG. (2.3) 


Naively, this volume is infinite since the integration extends all the way to the boundary at 
r = 0. Holographic renormalization proceeds in two steps. First, a regulator is introduced 
by extending the integration only to r = e. Second, counterterms are added to cancel e —> 0 
divergences, yielding a finite, physical answer. These counterterms must be built from data 
on the r = e slice to preserve covariance. This will include the induced metric on the r = e 
slice, 7 ij, and the Ricci scalar (and tensor) built from 7, i? 7 . This counterterm procedure is 
valid for any on-shell bulk action, which will in general be divergent due to integration down 
to r = 0, and applies to both the AdS d + ± volume and the free scalar action as we will review 
below. The counterterms for the volume renormalization are 


L 1 


L 2 



d(d-l)(d-2) 

ilog(e)v4#7 


7 


for d / 2 
for d = 2 


and for d = 4 only, an additional counterterm is needed: 

L 3 = - loge^iRijRV 


(2.4) 


(2.5) 


We will denote subtracted and renormalized quantities as V su b = V Ad s + / d d x[L\ + L 2 + L 3 ] 
and V ren = lim V sub . 

Hndices are raised and lowered using g( 0 ). 

2 This is related with the 011-shell value of the Einstein-Hilbert action for gravity with a negative cosmological 
constant as, Va<is = % ^°Ads ahM > where Sa<is = / d d+1 xVG(R + 2A). Notice that we use the conventions 

in [7], in particular the curvature conventions are such that AdS has positive curvature. 
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A free, massive scalar field <h(x,r) in this AdS background has action 


S = \f d d+1 xVG{G^d^d u ^ + M 2 <3? 2 ) (2.6) 

with equation of motion (EOM) the Klein-Gordon equation with the Laplacian associated 
with G^ v . Via the usual AdS/CFT dictionary, this scalar will be dual to some gauge-invariant 
CFT operator with dimension A given by M 2 = (A — d) A. A solution to the EOM resulting 
from the above action will, in general, have the form 

<h(x,r) = ?’ d-A 0(x,r) (2.7) 

= r d ~ A ((j) { o) + r</> (1) + ... + r 2A ~ d (l) (2A _ d) + r 2A ~ d log{r 2 )^ {2A _ d) ) + 0(r A+1 ). 


The logarithmic term appears only when 2 A — d is an integer, which we will assume to 
be the case, and reflects the existence of matter conformal anomalies in the dual CFT, i.e. 
conformal anomalies due to short distance singularities in correlators of composite gauge 
invariant operators [16]. All coefficients </>( n ) for n < 2A — d, as well as V’( 2 A-rf)) can be 
computed recursively by inserting this form of the solution into the EOM and expanding in 
powers of r. The result is that the leading coefficient, q i>( 0 ), and the coefficient 4>( 2A ~d) are 
not fixed by the EOM but all other coefficients are fixed by these two. The one coefficient 
we will need explicitly in terms of is the coefficient of the logarithmic term for the case 
2A - d = 2: 

0(2) = "Po0(O) + (| - 1)0(0) 1*0(2)]- (2-8) 

The action naively evaluated on such a solution is again infinite because of IR divergences, 
so holographic renormalization proceeds as before. The resulting counterterms are given by 


u = ^^-^(1 + 


1 


1 


£5 — 


d — A log e 


)d> 2 (x, e) 


2 ( 2 A -d- 2 ) 


V 7 ($(x,e)D 7 $(x,e) 


d -A 
2{d — 1) 


R[ 7]$0,e) 


(2.9) 


These counterterms are the only ones required when 2A — d < 2. The second term in L 4 , pro¬ 
portional to ] or l ( e \ , only appears in our examples when the scalar saturates the Breitenlohner- 
Freedman (BF) bound for AdS d + 1 , M 2 > —^ [17], for which 2 A — d = 0. In this case 
also L 4 alone suffices to renormalize the action. In the case that 2A — d = 2, the coefficient 
[2(2A — d — 2 )] _1 in L 5 is replaced with — | log(e) [7]. 

Renormalized correlators in the boundary theory can be computed from the renormalized 
action. According to the AdS/CFT dictionary, the leading coefficient of the asymptotic 
expansion will act as a source for the dual operator. Denoting this operator as O, the one- 
point function can be written 


(O) 


1 


SS r 


\/W) ^ 0 ( 0 ) ’ 


( 2 . 10 ) 
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which can in turn be written in terms of quantities on the r = e slice: 


(O) = lim 

e—>0 


/ J-_1 fiS sub \ 

\e A ^y6<f>(x,e)J 


( 2 . 11 ) 


For the cases we consider, if the scalar saturates the BF bound this will be modified: —► 

^ [17]. 

The structure of the counterterms and the renormalized one-point functions is most 
transparent in the radial Hamiltonian formalism, where the radius is playing the role of time 
[18, 19]. This is also the fastest way to arrive at the results reviewed above, since in this 
method covariance is manifest at all stages. In this formalism the central object is the radial 
canonical momentum, 


dL 

di’ 


( 2 . 12 ) 


where overdot denotes differentiation w.r.t. the radius p = — log?’. It can be easily shown to 
correspond to a regularized one-point function, 


< 213 > 

This expression is infinite as the regulator is removed. To extract the finite part in a covariant 
fashion we use the fact that the theory possesses a well-defined dilatation operator, 

fo = /A( 2 7ij H- + (A- rf )$T). (2.14) 

Let us expand the canonical momentum in eigenfunctions of 5d, 

7T = v^jd-A) H-b 7T(A) + 7T(A) ^g 7" 2 H-), (2.15) 

where the subscripts indicate dilatation weight of the terms, 

8d^(u) = ~ n7r (n)> £d7T( A ) = -A7T (A ), 8d^{ A) = ~^ 7T (A) ~ 2 ^(A)- ( 2 -16) 


Notice that the normalizable mode transforms anomalously. This expansion is a covariant 
analogue of the radial expansion (2.7). The explicit form of is found by inserting (2.15) in 
Hamilton’s equations and iteratively solving them by collecting terms with the same weight. 
This procedure determines all except 7r( A ) which is left undetermined, as in the discussion 
below (2.7). 

The renormalized 1-point function is now simply given by the term of weight A, as it 
should since the dual operator has dimension A, 


(O) = 7 T (A) . 

Furthermore, evaluating the regulating action on-shell one immediately gets 

S reg = - J d d x^n<S>, 


(2.17) 


(2.18) 
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and thus the counterterms are (minus) the terms involving the momentum eigenfunctions 
with weight less than A plus the logarithmic term involving tt\. This is a general feature 
of the counterterm action for any theory, namely the counterterms are always certain linear 
combinations of momentum eigenfunctions with weight less than A plus logarithmic countert¬ 
erms obtained by expanding the regulated action in eigenfunctions of the dilatation operator 
and using the basic relation (2.13) with 5 —> 5d, 

SnSreg = nS D $ + ^5Dlij (2.19) 

where we have reinstated the metric dependence. We refer to [18, 19] for the details and 
further discussion. 

3. Generalities for Probe Branes 
3.1 General Solution and Counterterms 

All the cases we consider in this section are SUSY probe D-branes embedded in AdS 5 x S 5 . In 
every case, the probe wraps some AdS m x S n with |m — n\ = 2 as required by supersymmetry 
[11]. We will focus on the cases with m—n = 2. In all these cases, a single scalar field describes 
the embedding. These probe branes are thus described by a scalar in a fixed gravitational 
background, only with an action more complicated than that of a free scalar. In general the 
non-linear Born-Infeld action would lead to new counterterms. We will show, however, that in 
the class of embeddings we analyze the counterterms are identical to that of a free scalar. To 
show this it is sufficient to show that the leading asymptotic structure of the field equations 
(i.e. the orders where the normalizable and non-normalizable modes appear) is the same as 
that of the free scalar, and establish that the on-shell actions agree up to order r m , the highest 
order at which divergences could appear in the action. From this follows (upon use of (2.13)) 
that the radial canonical momenta are the same up to tt(a) and thus the counterterm actions 
are the same. It also follows (and we will verify it explicitly) that the solution of the EOM 
has the same form as the free scalar solution to order r A . 

In three of the four cases we consider, the scalar is a coordinate on the S 5 transverse to 
the brane. The S n is a trivial cycle on the S' 5 , so the scalar, which gives the position of the 
S n in a transverse direction, can have a profile such that the S n “slips off’ the S 5 at some 
finite r. In the one exceptional case, for the D5 probe, the scalar is not a coordinate on the 
S 5 but is the position of the AdS 4, the brane wraps inside AdS$. This case is not included in 
the analysis of this section. 

The scalar “slipping mode” has mass-squared M 2 = — n (as we will show). Via the 
AdS/CFT dictionary the dimension of the dual operator is A = n. In particular, this means 
that d— A = (m—l)—n = 1, so in every case we consider the scalar’s asymptotic expansion will 
have a leading behavior of simply r. Recall also that the logarithmic term in the asymptotic 
expansion of the free scalar entered at order which here will be r n . The cases we consider 
are summarized in the following table of useful quantities: 
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Dp probe 

m = d + 1 

n = A 

2A — d 

cC z 

4 

D7 

5 

3 

2 

-4 

D5 

4 

2 

1 

9 

2 

D3 

3 

1 

0 

-1 


Given appropriate coordinates on AdS 5 , the metric of the AdS m subspace will have the 
same form as (2.1) but now /i, v = 0, ..., m — 1 and i,j = 0, ... ,m — 2. Let the S 5 metric be 

dQ 2 = d ( 3 ? 2 + sin 2 (<b)d0 2 + cos 2 (&)dfll (3.1) 

with dflg the standard metric on S 3 . The S n, s with n < 3 sit inside this S 3 . The DBI action 
is then 


S = J d m xVG cos n ($)Vl + Gt^d^dv® (3.2) 

where we neglect integration over the remaining n variables because the cases we are interested 
in have no dependence on these coordinates. The resulting EOM for $ is then 


0 = □<£ + n tan(<h) 


1 G^ v d^d v (G p<T dp<&d c $>) 

2 1 + G a Pd a <f>dp<!> 


(3.3) 


where □ is with respect to the AdS m metric. To order r A , this is identical to the EOM of 
a free scalar with M 2 = —n in a fixed AdS m background, which we explicitly show below. 
In the case of the free scalar the asymptotic solution has two undetermined coefficients, the 
source and the vev part, as we reviewed in section 2. This matches the fact that the field 
equation is a second order linear differential equation in r, so two pieces of boundary data are 
required in order to specify the solution. This is also consistent with holography: on the field 
theory side the form of the Lagrangian and the specific vacuum we consider (the vev’s) are 
the only information needed in order to specify the theory. We therefore expect that these 
facts will remain true for the asymptotic solution of the DBI field equation (3.3). Indeed, 
despite the higher derivative terms in the DBI action, the EOM (3.3) is still a second order 
differential equation in r (now non-linear) and we now verify this expectation explicitly. 

Consider a solution with the most general possible form: 

OO OO OO P 

<f>(x,r) = + ^*/j U) r J log{r) + ^2^2'I> {k: i ) r k log(r) 1 }. (3.4) 

i=0 j =0 k =0 1=2 


We are allowing higher powers of log(r), but inserting this solution into the EOM shows that 
generically truncating the series at any finite order p requires the t° a ll be zero. To see 

this, consider the leading term in the expansion of the EOM: 


0 = r[(m - 3)^(o) + 2T (0)2) + 21og(r)((m - 3)^( 0j2 ) +3T (0i3) )) + ... 
+ (P- 1) log(r ) p ~ 2 ((m — 3)T (0 p.!) +p’L (0jP) ) 

+ plog(r) p ~ 1 (m - 3)^ (0 ,p)] + Q(r 2 ) 
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Clearly, for any finite p. 4'( 0 , p ) = 0 implies 4'( 0) p-i) = 0 and so on so that all rlog(?’)^ do not 
appear when l > 2. A similar story happens for 4/( ljP p that is, all r 2 log(r) ; are absent for 
l > 2. Of course, an infinite number of higher powers of log(r) could appear (p = oo), but 
this would invalidate the assumption that the solution has a power series expansion. Special 
cases arise when the numerical coefficients multiplying the highest power of log(?’) are zero. 
For instance, m = 3 is a special case in (3.5) for which the coefficient of log(r) p_1 is zero 
(although in this case one still obtains ^(o.p) = 0)- These cases have to be analyzed separately 
and could result in higher powers of log(r) with coefficients determined in terms of the source 
and the vev coefficients. For m = 4,5 we verified that higher powers of log(r) do not appear 
up to order r 5 , that is, 4/( fc ^ = 0 for k = 2, 3,4, 5. When m = 3, 4b 22 ) and 4/ (2,3) are non-zero 
but are determined in terms of '(/’(o)- Notice that such higher order logs also appear in the 
Coulomb branch example in [17] and that both examples with higher powers of logs are cases 
where the BF bound is saturated. It is unclear to us whether this is a general feature (i.e. 
that higher order logs appear only when the BF bound is saturated) or it is specific to the 
examples we considered. 

What results from this discussion is that the form of the solution, up to the order we are 
interested in, is the same as that of a free scalar: 

A —1 A —1 

4*(x,r) = r[^2 4>{i)r l + ^ log(r)] + 0[r A+1 }. (3.6) 

i =0 j =0 

Inserting this into the EOM gives, to order r 2 , 

0 = (m — 3)V>(o ) r + ((m — 4)0^) + (m — 5 + (m — 4) log(?’))V’(i))r 2 + 0(r 3 ) (3.7) 

To this order, this is identical to a free scalar in AdS m with mass-squared M 2 = — n and 
d — A = 1. In particular, the leading coefficient, (j>(o)(x) is left undetermined by the EOM, as 
expected: this will be a source or will contribute to the vev of the dual operator, depending 
on whether the scalar saturates the BF bound. Notice also that for all m that we consider, 
V’(i) is zero. For m = 3,4, (3.7) determines the form of solution up to order r A . The value 
m = 5, however, requires the next order, r 3 . When m = 5, (3.7) requires '(/’( 0 ) = 0(1) = 0 in 
which case the r 3 term is 

0 = r 3 [D 0 (/>(o) + Tr(< 7 ( 2) )<^> (0) + ^{m - 5)^ 0) - 2(m - 5)0 (2) 

+ (7 - m)ip( 2 ) + 2 (5 - m) log(r)V’( 2 )] + 0(r 4 ). (3.8) 

This shows that </>( 2 ) is undetermined and ^( 2 ) is fixed by 0( O ) i R precisely the same way as 
the free scalar, (2.8), with d = m — 1 = 4. This is a nontrivial result: for all values of m we 
consider, the form of the solution is the same as that of a free scalar up to order r A . 

Since the solution is the same up to order r A so will the momentum dilatation eigenstates 
up to 7 T(a)- It follows that the divergences will be the same as that of a free scalar, and we 



explicitely verify this now. Expanding the argument of (3.2) we obtain: 

S = J d m xVG[ 1 + - n$ 2 ] 


n 


n 


1 


i 

+ (~ + ^)$ 4 - -( G ^ d ^ d ^) 2 - + 0 ( r 6 )]. 

1 z o o 4 

First notice that the leading term is simply the volume of AdS m , requiring the counterterms 

L\, L 2 and L 3 of the last section. 

The leading $ dependence is that of a free scalar in a fixed AdS m background with 
M 2 = —n, as advertised. As VG = and we care only about m < 5, all terms of 

order r 6 or higher cannot contribute to divergences. When m < 5 the <f> 4 terms, which go 
at leading order as r 4 , will not produce divergences, and hence in these cases the free scalar 
counterterms L 4 and L 5 are sufficient. When m = 5 these terms may produce logarithmic 
divergences, but these terms cancel against each other to order r 4 because n = 3: 

(“ + y)* 4 - \{G^d^d^) 2 - = [(~ + y) - l - J4)r 4 + 0(r 6 ) 

(3.10) 


(3.9) 


3.2 Supersymmetric Background Solution 


In the three cases with a scalar describing the slipping mode it is easy to check that with 
gij(x,r ) = r/ij an exact solution to the EOM preserving half the supersymmetries is 


$ = arcsin(cr) = cr + 'g' 7 ’ 3 + ■ • •, 


(3.11) 


3 

that is, a solution of the form (3.6) with = c, (j)^ = 0, f/v 2 ) = % and no logarithmic 
terms. This solution describes a D-brane extending up to 


T max — 


(3.12) 


c sets the mass of a fundamental matter field in the dual field theory. In the flat embedding 
space this embedding just describes a planar D-brane locate at a fixed distance away from 
the stack of D3-branes, 


X = -sin($) = c (3.13) 

r 

where X is one of the flat embedding space coordinates written in a spherical coordinate 
system. 

We will be mostly interested in calculating correlation functions, in particular the one- 
point function, in this particular background. In order to ensure that holographic renormal¬ 
ization gives finite answers for all correlation functions we need to make sure however to add 
all the counterterms for arbitrary boundary values of the fields, so that we get finite answers 
for expectation values in the presence of sources. While the counterterms and the background 
solution are very similar in all the cases, the physical interpretation as well as the details of 
the procedure are quite different in all the cases and hence we will turn to the examples one 
by one. 
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4. The D3/D7 System 

4.1 The System 

In the supersymmetric embedding of a probe D7 in AdS§ x S' 5 , the D7 fills the 5D spacetime 
and wraps an S 3 inside the internal space, that is, the probe wraps an AdS 5 x S 3 and hence 
m = 5, n = 3, and the scalar’s mass is above the BF bound. This system was proposed 
in [9] as a way to introduce flavor, that is fundamental representations, into the AdS/CFT 
correspondence. The particular supersymmetric case we are studying is dual to the M = 4 
SYM theory coupled to a fundamental hypermultiplet that preserves N = 2 supersymmetry. 
Already in [9] it was noted that the gravity side has a non-zero subleading term ^( 2 ) in 
its asymptotic expansion, even though from the field theory side it is clear that no vev is 
allowed. Some arguments were given that on the gravity side the vev really is zero. Starting 
from [ 20 ] it has become common practice to read off the vev from the subleading term in the 
flat embedding space coordinate X = isin(<h) which on the background solution is constant 
and indeed has no subleading term. The absence of any rigorous derivation of the correct 
procedure to calculate the vev became apparent with the work of [ 21 ] where different answers 
were found depending on what coordinate system was used. In their work a very natural 
choice was again “apparent” that yielded the expected answers, but the question arose how 
one would find the right coordinate system in less symmetric cases. Here we will apply 
holographic renormalization to that problem and will find an unambiguous answer. 

4.2 Counterterms 

In addition to the counterterms we derived for general m and n above, only in the case of 
the D7 brane we face an interesting new subtlety: the possibility of finite counterterms. A 
counterterm of the form 

Lf = oiy/' 7 <h(x,e ) 4 (4.1) 

will not introduce divergences for any a , but will change the on-shell action by a finite 
amount depending on the free coefficient a. Different values of a correspond to different 
renormalization schemes and some correlation functions will contain scheme-dependent terms. 
As first pointed out in [17] in the special case that one is interested in correlation functions in a 
supersymmetric background, one can fix the finite counterterms by picking the unique scheme 
in which supersymmetry is preserved. In the supersymmetric scheme not just the divergent 
pieces of the on-shell action but the action as a whole has to vanish when calculated on the 
supersymmetric background solution. Otherwise the ground state energy in the dual field 
theory would not be zero and hence supersymmetry would be broken. 

For the D7 case the counterterms read 

Li = ~\Vl (4.2) 

U = " 48 ^^ 
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L 3 = - log eyfy^RijR' 3 - ^1? 7 ) 

^4 = -y/i$ 2 {x,e) 

L 5 = -^log(e)^$(x,e)(p 7 + ^R 7 )4>(x,e) 

For the supersymmetric embedding, gij(x,r ) = rjij and the curvature on the slice is 
zero so only L\ and L 4 contribute. To fix the finite counterterm we only need to plug the 
background solution (3.11) into the action. On the background solution the action exactly 
reduces to 



Indeed the action picks up a non-vanishing contribution of — 7 -^-b = +x from the IR 

boundary. The UV divergence is ensured to cancel by the counterterms we introduced above. 

4 

L 4 gives an additional finite piece, L 4 = ... + All in all we see that in order to set 

the on-shell action to zero we need 


a = — 


_ 5 _ 

12 ' 


(4.4) 


4.3 One-point function 

Holographic renormalization tells us that the vev is obtained as the e —> 0 limit of 

1 SS sub 


( 0 ) = 


e 3 ^/7 54>(e) 


In this case, 


5S reg 5L 


<5<h 




= —3 v /5cos 3 (4>(e)) 


4>'(e) 


V 7 ! + e 2 4 >'(e) 2 


(4.5) 


(4.6) 


This then contributes divergent and finite parts to the vev 


(O)reg ~ -4jr “ 30( 2 ) l°g( e ) + (20( O ) ~ ~ '0(2)) + 0(e 2 ) (4.7) 

Adding in addition the variations of the counterterms we finally obtain for the vev 

(O) = log(e)(— 2 ^( 2 ) — -0(o)-Ro — ^ 00 ( 0 )) + (3 0(o) — 20( 2 ) — 0 ( 2 )) + 0(e 2 ). (4.8) 

The coefficient of the log(e) is zero via the EOM, (2.8), and inserting the value of 0( 2 ) from 
( 2 . 8 ) into the finite piece then gives 


(O) 


— 20(2) + 



+ -j7^0(O) + n O0(O)- 


(4.9) 
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The last two terms come from the curvature of the induced metric on the r = e slice and can 
in fact be eliminated by adding a finite counterterm proportional to the matter conformal 

3 

anomaly. On the background solution sin(<E>) = cr we have = c and <M) = so indeed 

(O) background — O' (4.10) 

The finite counterterms required by supersymmetry set the vev to zero despite the appearance 
of a non-vanishing subleading term in the asymptotic expansion of the scalar field. 

5. The D3/D5 System 

5.1 The System 

In the supersymmetric embedding of a probe D5 in AdS 5 x S' 5 , the D5 fills an AdS 4 subspace 
of the 5D spacetime and wraps an S 2 inside the internal space, that is, the probe wraps an 
AdS 4 x S 2 and hence m = 4, n = 2 and the scalar’s mass is above the BF bound. This 
system was introduced in [10] and analyzed in detail in [22], The dual field theory in this 
case is a defect conformal field theory. Again the probe brane in the bulk corresponds to the 
introduction of fundamental matter in the field theory, but this time the matter only lives on a 
3d defect in the 4d field theory. While the defect breaks translation invariance in the directions 
orthogonal to the defect, the field theory remains conformal as long as the matter is massless. 
The subgroup of the full conformal group that is preserved is the SO( 3,2) that leaves the 
position of the defect invariant. This symmetry in the bulk demands the worldvolume to be 
AdS 4 . The supersymmetric background solution once more describes adding a mass term to 
the fundamental matter. A new possibility that arises in this case is to turn on the scalar that 
corresponds to the embedding of the AdS 4 inside the AdSs [11]. This deformation breaks all 
supersymmetries. On the field theory side it was argued in [11] that one turned on a vacuum 
expectation value for the defect field, as we will verify explicitly below. We will also see that 
the energy of the configuration in non-zero. On the field theory side one certainly expects 
such configurations to eventually relax to the supersymmetric ground state. On the gravity 
side, this process should be mapped to some instability of the brane embedding that drives 
the vev to zero. It would be interesting to analyze the stability of this embedding. This 
would provide an example of an unstable (or perhaps metastable) configurations that exists 
and can be studied both at weak and strong coupling. 

5.2 Mass deformation 

In order to cancel the UV divergences we need to add the following counterterms: 

L \ = ( 5 - 1 ) 

L 2 = 

L 4 = -y/i$ 2 (x,e) 
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L 5 = -^v'T : $(x,e)(P 1 + ^R 1 )$(x,e) 

No option to add finite counterterms involving only and arises. Indeed with this set of 
counterterms the on-shell action on the background solution is zero already, 


Sreg — 



(1 — c 2 r 2 ) 3 / 2 
3 r 3 



(5.2) 


The IR term, with r max = —, is identically zero. 

With these counterterms it is again straightforward to calculate the vev. Setting to zero 
the terms involving i ? 7 and D 7 one simply obtains 


(O) = -</>(!)• 


(5.3) 


Again, the supersymmetric background solution which has 4>{i) = 0 has no vev, as expected. 


5.3 Vev deformation 


In this subsection we study the probe RG flow obtained by turning on a scalar that controls 
the embedding of the AdS 4 part inside the AdSs while keeping the brane wrapped on the 
maximal sphere inside S' 5 , T = 0 [11]. For simplicity in this case we will only calculate the 
vev on the background solution; no problems are anticipated in once more working out the 
full holographic RG program for arbitrary correlation functions. Looking for an embedding 
X 3 = x(r) the induced metric on the 4d part becomes 

ds 2 = -^2 (— dt 2 + dx\ + dx 2 + (1 + (x') 2 )dr 2 ) (5.4) 

and the corresponding action is 


s = /^y ttw. 


(5.5) 


The most general solution to the equations of motion with the ansatz that X 3 only depends 
on r is 


cr 


x = 






so that asymptotically 
This brane extends up to 


x = ?’o H—r° + ... . 


1 


m.n.T. 


a/4 * 


(5.6) 

(5.7) 

(5.8) 


The solution with c = 0 is the supersymmetric AdS 4 x S 2 defect. What happens at r = r max 
for the vev deformation is quite distinct from what we have encountered in the other exam¬ 
ples. While before a sphere was shrinking in the internal space so that the brane smoothly 
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terminates, for the vev deformation the internal sphere has constant size. To understand the 
physics of this embedding it is useful to note that close to r max we simply get 


/ WV m ax 

X = 


2 v/2 y/fmax T" 


X = 


-Vr 


— r + 


(5.9) 


Instead of ending at r = r max the string turns around and hits the boundary again. The 
probe brane no longer is dual to the defect conformal field theory of a single D5 brane defect, 
but with a D5 and anti-D5 defect seperated by a distance d = 2 
The on-shell action can be evaluated exactly 


Sreg — 






3/4 v / 7rr(13/8) 
15r(9/8) 



(5.10) 


where the finite term comes from the IR boundary, 
action the counterterm 


It follows that we should add to the 

(5.11) 


With this counterterms the action is not just finite but zero for the supersymmetric embed¬ 
ding, and is negative for the non-super symmetric embedding. For the general x(r) solution 
we get for the vev 


in\ = 1 1 SSsub = 1 x ' 

e 3 y/j 6x e 4 yj\ + (x'Y 


(5.12) 


This agrees with the expectation that a non-vanishing c corresponds to a vev-induced flow 
once one realizes that ^ is dual to an operator of dimension four [22]. The negative energy 
does not signal an instability of the original D5 brane defect since, as we noted above, the 
vev deformation really is dual to a system with 2 defects, one due to a D5 and one due to 
an anti-D5. It is not too surprising that latter system can lower its energy by moving out 
onto the Higgs branch. Since the constant c is completely fixed by the distance d between the 
defects we expect this configuration to correspond to the stable ground state of the system. 


6. The D3/D3 System 
6.1 The System 

The D3/D3 system is another example that can be treated with the same methods. In this 
case, the probe D3’s wrap an AdS$ x S 1 inside AdS$ x S 5 , m = 3 and n = 1. This system was 
analyzed in [12]. It corresponds once more to a defect conformal field theory, this time with 
the fundamental hypermultiplet localized on a 2d defect in the 4d field theory. The interesting 
new aspect of this system is that the mass of the slipping mode saturates the BF bound. This 
gives rise to a slightly different structure of counterterms in holographic renormalization; it 
also leads to an interesting reversal in the role of the vev versus mass as we will see below. 

From the boundary theory point of view it is also clear that this system is different. In all 
the other cases we had a higher dimensional brane intersecting the D3 brane, so that in the 
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field theory limit the worldvolume fields on the higher dimensional brane decoupled together 
with gravity. The only dynamical fields remaining were those on the D3 branes and on the 
intersection from the 3-7 or 3-5 strings respectively. For the D3/D3 system the M = 4 SYM 
on both D3 stacks remains dynamical, so the field theory is really the theory of two J\f = 4 
theories coupled to each other via a bifundamental hypermultiplet living on a 2d defect. In 
particular, the spacetime separation between the two stacks of D3 branes can be viewed as 
a vev for one of the additional D3’s worldvolume fields and should no longer be interpreted 
as a parameter (as in the D5 or D7) case, but as a vev of a dynamical held. We will see 
that holographic renormalization correctly accounts for this effect. It is precisely separating 
the two stacks in spacetime that is described by the supersymmetric background solution we 
studied for general m and n, and we will see that unlike in the other cases for m = 3, n = 1 it 
is a vev instead of a mass deformation. Implementing holographic renormalization properly 
we also calculate the two-point function in this case, yielding a different answer from the one 
originally advertised in [12]. 

6.2 Counterterms and One-Point Function 

The interesting change in roles between source and vev can already be observed from the 
form of the asymptotic solution which in this case reads 

$ ='K^(o) + ^(o) log(j) + • • •) (6-1) 

for which the first term is normalizable and will set the vev, while the second term is non- 
normalizable and will correspond to a source in the dual theory. On our supersymmetric 
background solution the source term is zero; we will calculate the vev momentarily. For 
m = 3, n = 1 the counterterms become 

L i = -\\/l ( 6 - 2 ) 

L 2 = ^ log ey'yiiy (6.3) 

l 4 = ^v/7^ 2 (a;,e)(l + j^) (6.4) 

No finite counterterms can arise and it is again easy to verify that the on-shell action for our 
supersymmetric background vanishes with these counterterms already. 

In this case, the vev is the e —> 0 limit of 

1 SS su b 

v/7 Mo¬ 
using the counterterms it is again straightforward to obtain 

(O) = 4> {0) ( 6 . 6 ) 

where we again set to zero the terms involving the curvature on the slice. The vev is thus 
indeed determined by the coefficient of the normalizable mode; on our supersymmetric back¬ 
ground we simply have (O) = c. As expected, the brane separation shows up as a vev and no 
longer as a mass term. 


(6.5) 


(0) = 


log(e) 
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6.3 Two-point functions 

We discuss in this section the computation of correlation functions, n-point functions can be 
obtained from the exact one-point function by further differentiating (n — 1) times w.r.t. the 
source. It follows that in order to obtain the two-point function it is sufficient to solve the 
linearized equation of motion. 

The simplest case to consider is the two-point function for the theory specified by the 
embedding $ = 0. We first discuss the 2-point function of the operator dual to the scalar field 
<h. This can be obtained by linearizing the field equation (3.3) around <J> = 0, and Fourier 
transforming in the spatial directions we get 

r 2 <F" _ r & _ (fc2 r 2 _ = o ( 6 . 7 ) 

This is equal (as it should be) to the linearization of the field equations of a scalar of mass 
M 2 = — 1. The solution of this equation that is regular in the interior is 

&(r,k) = rKo(kr). (6.8) 


Expanding for small r we get 


®(r,k) = Vfo )(k)r 


1 k 2 

log r + - log — 

Z jjj 


(6.9) 


where ^(o) (&) represents the overall normalization of the solution, and fi 2 = 4e 27 (7 is the 
Euler constant). Proper units are restored using the AdS radius. It follows that 

Vfo) (*) = Vfo) (*) l log ( 6 - 10 ) 

Z fZ 

{0(k)0(-k)) = -fp- = -Uop- 2 (6.11) 

Fourier transforming we get the renormalized version of 1/x 2 (see appendix) 


(O(x)O(O)) = A (P) (6.12) 

(with m = 1), as is appropriate for the two-point function for a scalar operator of dimension 
one. 

We now discuss the computation of 2-point function which are dual to fields parametrizing 
the flunctuations of AdS% in Ad,S§. The mass of the fluctuations is [12], 


m 2 = (l — 1)(Z — 3) 


(6.13) 


with l a positive integer. One should distiguish between the fluctuations with mass 


2 « 

m 2 > —- + 1 = 0 


l > 3 


(6.14) 
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in which case there is a single branch of dual operators with dimension, 

A = / — 1, l> 3 (6.15) 

and fluctuations with mass 

d 2 d 2 

= -1 < m 2 < +1 = 0, =4- 1<£<3 (6.16) 

In this case there are two inequivalent quantizations and we have two branches of operators. 
The l = 2 correspond to m 2 = —1, so the fluctuations saturate the BF bound. We discussed 
the computation of the 2-point functions for this case in the beginning of this section. The 
l = 1 and l = 3 cases correspond to massless fluctuations. The l = 1 case is special: the 
corresponding boundary operator has dimension zero and thus saturates the unitarity bound. 
From the bulk point of view this case utilizes the A_ branch and the roles of source and vev 
are exchanged. The computation of correlation functions for such cases has been discussed 
in [25]. When the operators saturate the unitarity bound, however, the normalization of the 
(holographically computed) 2 -point function turns out to vanish (set A = d/2 — 1 in (2.19) 
or (2.21) of [25]) even though the corresponding QFT 2-point function is generically non-zero 
indicating subtleties with the normalization of these operators. Leaving this special case aside 
we now discuss the computation of 2-point functions for l > 3. 

In our analysis so far we have constructed the counterterms required to render the on- 
shell action finite and then derived the 1-point function. This can also be done for the case at 
hand and this would uniquely fix (up to finite scheme dependent terms) all boundary terms. 
We should also note that another (perhaps more fundamental) requirement on the action is 
that it is stationary when the field equations hold and appropriate boundary conditions are 
imposed. This also fixes all boundary terms to be the ones determined by finiteness [23]. 

In this section we are only interested in computing 2-point functions, however, so instead 
of first computing the counterterms that would render finite any correlation function we 
will follow the simpler route developed in [19]. 2-point functions are determined by solving 
the linearized fluctuation equations and to renormalize we need to perform a near-boundary 
analysis of the fluctuation equations. 

The fluctuation equations can be obtained from (3.16) of [12] 3 . In our coordinates these 

read 

r 2 w'( — rw[ — ( k 2 r 2 + m 2 )wi = 0 (6-17) 

where we Fourier transformed along the boundary directions. In the recent work [24] the 
fluctuation equations were analyzed for D3 branes separated by a distance L and shown to 
be of the hyper geometric type. One can thus readily derive the 2-point functions in this more 
general setting but for simplicity we restrict to the L = 0 case. The fluctuation equation in 
(6.17) can also be obtained from (4.10) of [24] in the limit L —■> 0. The solution of (6.17) that 
is regular in the interior is 

wi(r,k) = rK {l _ 2 ) (kr) (6.18) 

3 We employ the w notation of [12], whereas [24] uses w. The two are related by w = rw 
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We now discuss how to extract renormalized 2-point functions from the solution of the 
fluctuation equation. We will follow the Hamiltonian formalism described in the introduction. 
It is useful to introduce a new radial coordinate p = — logr. In these coordinates the AdS^ 
metric reads 

d-s 2 = dp 2 + e 2p dx l dx l (6.19) 

and the canonical momenta are given by 4 

7T/ = Wl (6.20) 


where overdot denotes differentiation w.r.t. p. The fluctuation equation (6.17) becomes 


7 Ti + 27T; — (£ + m 2 )wi = 0 


( 6 . 21 ) 


where we defined £ = e~ 2p k 2 , which comes from the Fourier transform of at the regulating 
surface r = const and 7 ^ = e 2p 5ij is the corresponding induced metric. 

Since we are interested in 2-point functions, we need to know the radial canonical mo¬ 
mentum as a function of the induced field to linear order. Covariance fixes their relation to 
be the form 

ti 7 = MOW- ( 6 - 22 ) 

The regularized 2-point function (with the regulator being a small constant value of the radial 
coordinate r = e) is then given by 


<0,_i(A;)0,_i(-A;)> e = 


1 

g2(Z—2) 


M) 


1 

g2(Z—2) 


(^(l - 3) + ke 


K(i-2)(ke)J 


(6.23) 


In the limit e —> 0 these correlators diverge. Our task is to extract correctly the finite part. 

As reviewed in section 2 , the renormalized 1 -point function in the presence of sources is 
given by the part of the canonical momentum with dilatation weight equal to the dimension 
of the dual operator [18, 19]. To extract this part we need the asymptotic form of /)(£). To 
this end, we insert ( 6 . 22 ) in ( 6 . 21 ) to obtain 


fl + fi + 2 fi — (£ + m 2 ) — 0 


(6.24) 


From the general discussion in the introduction we know that the //(£) has an expansion of 
the form, 

fliO = fi( 0) + fi( 2 ) 4 -F f 1 ( 21 - 4 ) + f 1 ( 21 - 4 ) log ?’ 2 -I- (6.25) 

Furthermore, in our case the dilatation operator in (2.14) is just equal to the radial derivative, 


$D = d P 


(6.26) 


4 Strictly speaking, the momenta are densities and one should include a factor of as in the formulae in 
the introduction. The equations are simpler without these factors, however, so we will work with (6.20). 
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This follows from the fact that the background is AdS 3 (insert 7 ^ = e 2p 5ij, <f> = 0 in (2.14)). 

Inserting the expansion (6.25) into (6.24) and organizing the terms according to their 
dimension one obtains all coefficients except for f 1 ( 21 - 4 ) which is left undetermined. Let us 
explain this computation. In (6.24) the mass term has dimension 0 and £ has dimension 2, 
so fi( 0 ) and /^ 2 ) are different from the remaining fm p )- We get 

1 + 3 : f m = 1 - 3 , f l{2) = (6.27) 

1 = 3 : / 3(2 ) undetermined, / 3 ( 2 ) = 

For the remaining coefficients one gets a recursive relation, 

1 P_1 

fl(2p) = 2 ) fl(2(p-n))fl(2n) (6.28) 

' n=l 

^ 1-3 

f 1(21-4) undetermined, f 1(21-4) = - ^ fi(2(l-2-n))fl(2n) 

n =1 

Having determined the asymptotic expansion we now extract the renormalized 2-point 
function as 


P¥=l~ 2 : 

p = l — 2 : 


(0;_i(£:)0 ; _i(-&)> = - lim f (/ - 3) + ke ^ f l{2p ) - f 1 ( 21 - 4 ) loge 2 


K(i- 2 )(ke 


p =0 


Let us work out explicitly the first few cases (/U 2 = 4e 27 as before): 


(6.29) 


1 = 3: (0 2 (k)0 2 (-k)) = h 2 log^ (6.30) 

2 /i z 

Z = 4: (0 3 (fc)0 3 (-fe)) = -i(fc 2 ) 2 log^ 

8 

l = 5 : <O 4 (A0O 4 (-A0> = y^^Oog ^ 


which is the correct behavior for two-point operator of dimension 2, 3 and 4, respectively. 
(The last term in the r.h.s. of the l = 5 case is scheme dependent. It contributes a contact 
term in x space). Fourier transforming to x space (using the results of the appendix) we find 


(Oa(x)Oa(O)) 


2vT(± + v) f 1 \ 
Tt d / 2 T(v) V'C 2A /^J 


(6.31) 


where the subscript R indicates that the expression is renormalized and v = A — d/2. The 
normalization of the 2-point functions coincides with the one determined in [26] (as it should). 
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By construction, the subtractions in (6.29) can be implemented by means of counterterms. 
The explicit form follows from (2.18) and is given by 

= d 2 x^wi I Y fl(2p) + log e 2 /z(2Z-4) ) W* + c.c. 

4 l >3 J \p =0 ) 

= \ Y f d 2 Xy/l ({l - 3 )wiwt - 2 ( /_ 3) ^D>r + • • + C.c., (6.32) 

(when l = 3 there is only the logarithmic counterterm). The renormalized action is given by 
the sum of the bulk action and this boundary term. The action derived in [12] should be 
viewed as a “bare” action. Adding to it appropriate boundary terms leads to a renormalized 
action and the corresponding renormalized correlators exhibit power law behavior. 

The fields discussed here are dual to supersymmetric operators that saturate a BPS 
bound [12]. A realization of these operators in terms of defect fields has been proposed in 
[12] and they contain massless scalar fields. As discussed there this implies that in general 
there will be strong infrared effects. The holographic computation discussed here and the 
fact that these operators are supersymmetric suggest that appropriate IR renormalization 
would lead to 2-point functions exhibiting power law behavior at weak coupling. We leave 
the investigation of this issue for future work. 


Appendix: Fourier Transform 

In differential regularization [27], 1/x 2 is represented by 


= -□log 2 .x 2 m 2 


(6.33) 


where the subscript R indicates that this is a renormalized expression and w? is the renormal¬ 
ization scale. The renormalized expression differs from 1/x 2 by the infinite term 5(x) log x 2 m 2 
localized at x = 0. 

We now compute the Fourier transform following [28]. Integrating by parts once and 
after some manipulations one arrives at 


XI „„ikx 1 


,.2 „2^,2 _ • 4* f i2„Akx X ^ 


I(k) = / d z xe ikx -B\og 2 x 2 m 2 = -i-£ / d 2 xe lkx ^ log x 2 m 2 
J 8 2 J x z 

= -4 ((?) F(a) ) L=o 


(6.34) 


where u 2 = 4e 27 m 2 and 


o2a+l r 

F(a ) = —a-e _2 “ 7 / rdrdO- 

2tt J 


rp2a+2 


(6.35) 
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The integral can be readily computed, 


F(a) 


= -27a r ( X ~ a ) 

' r(i + o) 


l + 0[a] 3 


It follows 


I(k) = — 7 r log 


H 2 


(6.36) 


(6.37) 
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